General Relativity and Gravitation manuscript No. 

(will be inserted by the editor) 



o 

(N 



O 

U 



O 
O 



O 

o 



The Relevance of the Cosmological Constant for Lensing 
Mustapha Ishak • Wolfgang Rindler 



Received; 26 November 2009 / Accepted: 14 March 2010 



Abstract This review surveys some recent developments concerning the effect 
of the cosmological constant on the bending of light by a spherical mass in 
Kottler (Schwarzchild-de Sitter) spacetime. Some proposals of how such an 
effect may be put into a setting of gravitational lensing in cosmology are also 
(~| , discussed. The picture that emerges from this review is that it seems fair to 

Qh' assert that the contribution of A to the bending of light has by now been 

Q , well established, while putting the A light-bending terms into a cosmological 

i—^ • context is still subject to some interpretation and requires further work and 

C/3 , clarification. 



Keywords bending of light • lensing • cosmological constant 



> , 

■"nJ" I 1 Introduction 



In view of the present candidacy of the cosmological constant yl as a possi- 
ble cause of the observed acceleration of the universe, any non-cosmological 
manifestation of A could be of interest. Already in his 1923 treatise The Math- 
ematical Theory of Relativity [Sec. 45 in [1]]. A. S. Eddington examined one 
possible such manifestation, namely the contribution of A to the perihelion 
advance of the planets. Soon thereafter he applied his results to the special 
case of Mercury, and found that if A were greater than 5 x 10~''^cm~^ it would 
■ increase the centennial advance of the perihelion from 43" to more than 44", 

H I which would have been detected by the astronomers, see also [2]. Eddington 

never seems to have made a similar calculation for the other then-crucial effect 
of GR, the bending of light. Possibly in view of Einstein's later renunciation 



Department of Physics, The University of Texas at Dallas, Richardson, TX 75083, USA 

Tel.: 1-972-883-2815 

Fax: 1-972-883-2848 

E-mail: mishak@utdallas.edu 

E-mail: rindler@utdallas.edu 



2 



of the yl-term, interest in this matter faded. In any case, when compared to 
the cosmologically determined value of about 10~^^cm~^, Eddington's upper 
bound strongly suggests that a "locally" detectable yl-effect could only be 
expected on cosmologically significant length scales. 

In modern times the problem was taken up again apparently for the first 
time in a much-quoted 1983 paper by N. J. Islam [3]. There, the author re- 
derived Eddington's A limit of 10~^^cm~^ from the perihelion advance of 
Mercury; but, more importantly, he comes to the conclusion that A has no 
influence on the bending of light. The case for this noninfluence of A was then 
remade and reaffirmed by various other authors. N. J. Islam's argument was 
very simple: In the Kottler metric [4], which is the modified version of the 
Schwarzschild metric when A is included in the field equations, A modifies the 
orbits of all massive particles, but drops out of the coordinate orbital equation 
for photons. 

However, in a 2007 paper by the present authors [5] it was pointed out 

that the coordinate differential equation for the orbit and its integral are only 
half the story. The other half is the metric itself. It is only the metric that 
determines the geometric and measurable properties of the coordinate orbital 
equations. When that is taken into account, a quite different picture emerges: 
A does indeed contribute to the bending of light. 

Apart from its fundamental aspect, this matter is of potential interest to 
the many people involved in gravitational Icnsing where the approach to the 
problem is generally somewhat different from the basic geometrical one of [5] . 
Nevertheless, most authors seem to have accepted the existence of Tl-bending in 
principle, though some have expressed reservations, and others have obtained, 
by different methods, somewhat different formulae as we will discuss. 

2 Light bending in Kottler (Schwarzschild-de Sitter) spacetime 

2.1 The contribution of the cosmological constant to the bending of light 

As is well known, the Schwarzschild metric is the unique solution (up to co- 
ordinate transformations) of Einstein's original vacuum field equations, when 
applied to the vacuum between a spherically symmetric central massive object 
and a possibly surrounding concentric mass distribution (or simply an infinite 
vacuum). By Birkhoff's theorem [6], any purely radial and spherically sym- 
metric motion of the central object or the surrounding masses has no effect on 
the staticity of the intervening vacuum. 

When Einstein's equations are extended by the so-called cosmological or 
yl-term, the unique successor to the Schwarzschild metric was first obtained 
by Kottler [4]: 



(1) 



where 



2m Ar^ 



(2) 
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Fig. 1 Schwarzschild and Schwarzschild-de Sitter geometries. 17^ is the Flamm paraboloid 
representation of a central coordinate plane in Schwarzschild; 17^ is the corresponding surface 
in Schwarzschild-de Sitter; is an auxiliary plane with an r,<f> graph, , of the orbit 
equation (6). The curves and are the vertical projections of onto S'^ and E'^ , and 
represent the true spatial curvature of the orbits. 



and where units are chosen so as to make c = G = 1 , as wiU be done throughout 
this paper unless otherwise stated; m is then the mass of the central object 
expressed in length units, the only units remaining. When A = Q, Kottler's 
metric reduces to Schwarzschild's, and when m = it reduces to the static 
version of de Sitter space. These limits are also approximated by the Kottler 
metric when r is small and when r is large, respectively. By an extension of 
Birkhoff 's theorem [6] , any purely radial and spherically symmetric motion of 
the masses involved leaves the staticity of the Kottler metric unchanged. 

The geometric modification that Kottler brought to Schwarzschild is highly 
significant for our purposes. Because of the spherical symmetry of the spatial 
part of both these metrics, it is clear that any orbit of a particle or a photon 
must lie entirely in one of the central coordinate "planes" typified by the 
equatorial plane 9 = 7r/2. In fact all orbit calculations are usually made in 
this representative plane. Of course, although the locus — it/ 2 is plane in 
the coordinate sense, it is not plane when measured with rulers. The ruler 
geometry of this locus is determined by the spatial part of the metric (1), 
which for the Schwarzschild case {A = 0) is 




(3) 



and for the dc Sitter case (m 



0) is 



(4) 
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As is well known, (3) has the intrinsic geometry of the upper half of the so- 
called Flamm Paraboloid [7] which, when displayed as a surface of revolution 
in ordinary Euclidean space, has the equation 



2 



z 



8m{r-2m), (5) 



z being the distance along the axis and r the distance from the axis. The 
corresponding intrinsic geometry of the metric (4) is that of the lower half of a 
2-sphere of radius y/^/A, with the origin at the southpole, (j) being longitude 
and r again the distance from the axis. In the Schwarzschild case the range 
of r (region of staticity) is bounded below by r = 2m, the waist circle of the 
paraboloid. In the de Sitter case it is bounded above by r = -^/S/yl, the equator 
of the sphere. 

The spatial geometry of the Kottler metric (1), being essentially a combi- 
nation of (3) and (4), is illustrated in Fig.l. For small r it approximates the 
Flamm funnel down to r « 2m, while for large r it widens into an approximate 
de Sitter sphere up to r « ^/sJA. The precise r- values of these boundaries are 
given by the two positive zeros of f{r). These correspond to event horizons, 
but will not concern us here, since the region of interest lies well in between. 

The work of [5] is based on the well-known first approximation orbital 
equation for photons in the 9 = tt/2 plane of the Schwarzchild metric 

1 sinid)) 3m / 1 , A 

which, as mentioned earlier, applies equally in the Kottler metric. The orbits 
are parameterized by i?, the distance of the zeroth-order (straight-line) approx- 
imative solution from the center (see Fig. 2). The relation of R to the physically 
meaningful coordinate distance (area distance) rg of closest approach is ob- 
tained by setting (/> = 7r/2 in (6): 

Now, in Schwarzschild space it is easy to obtain the total deflection angle suf- 
fered by a ray of light coming from infinity and going to infinity by letting 
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r — > oo in (6). At infinity, the coordinate plane = 7r/2 is essentially Eu- 
clidean, and the coordinate angle </> is also the measured angle (see Fig. 2). 
Not so in Kottler space. There r — > oo is impossible because of the spherical 
nature of the geometry (sec Fig.l). The approach adopted in [5] is to find the 
angle if' measured by an observer at an arbitrary point {r,4>), between a ray 
coming radially from the center and the ray skirting the lens (see Fig.2). This 
angle is given by the standard formula 

for the angle between and (5*. It is precisely here where A comes into play, 
via the metric coefficients gij of the Kottler metric (1). 

For an r-value sufficiently far from the lens for (p and tp to be small (see 
Figure 2), the results of [5] (especially Eqs. (13) and (16)) can be shown (by 
a sequence of Taylor approximations, working to first order in (j), ip, m/R and 
Ar'^) to be equivalen to 

, 2m AR^ 
'^ = '^+^- 6{2m + cpRy 

This would seem to definitely indicate a contribution of A to the bending of 
light. [The term in the corresponding Eq. (17) of Ref. [5] is superfluous 
and innacurate since the expansion was not done to the second order but see 
equation (18) below] As can be seen from Fig.2, the one-sided actual deflection 
angle is ijj — 4>. Thus, in the Schwarzschild case where r — > oo makes sense, 
r — > oo implies ip — > and then (9) gives the familiar Einstein deflection. 

A particularly transparent application of Eq.(9) is given in [9]. It concerns 
an Einstein ring, a well-known example of gravitational lensing, where a con- 
centrated source directly behind the lens is seen by an observer directly in 
front of the lens apparently spread into a ring. The case most easily dealt with 
by formula (9) is when the source and observer are equidistant from the lens, 
so that in Fig.2 the observer is at i/i = and consequently, by (9), 

, 2m AR^ 

ib = . (10) 

^ E 12m ^ ' 

This, in principle, appears to offer a method of measuring the yl-contribution 
if we know m and ro, the radius of the lens - which is related to R by Eq. (7) 

Of course, lensing on the cosmological scale is not done by an observer 
at rest relative to the source. But any radial velocity V of the observer would 
simply increase the entire ^p of Eqs. (9) and (10) by the usual aberration factor 
1/(1 + V)/{1 — V) of special relativity, and thus leave the relative contribution 
of A unchanged. Neverthless, in a universe with cosmological constant yl, the 
velocity V is yl-dependent, and the cosmological expansion effect of A (via 
aberration) can exceed its above discussed geometric effect by many orders 
or magnitude. For example, in a simple de Sitter universe, where V « Hr = 
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A/Zr, the abberation factor is ^ [1 + ^J A/2)r) and then the dominant yl-term 
in (10) becomes 

2mr [a ^ [a PIrwrA 

\ — = ii\ — = \ , (11) 

ii V 3 V 3 V 3 ' ^ ' 

where we have used Eq. (6) with = to relate r, R, and m. Since Va IS very 
much bigger. Since A is small, root A is much bigger, and so the aberration 
effect of A overshadows its geometric effect. 



2.2 Higher order mass terms 

ft is of interest to compare the yl-terms to higher-order mass terms, for which 
purpose we must extend the calculation above to higher orders. Proceeding in 
the usual way, see for example [10,11], one writes 

u = uo[sin{(f)) + {muo)Sui + {muo)'^Su2\ (12) 

where u = - and uq = ^. Substituting this into the null orbit ODE, e.g. [8, 
12] 

— +u = 3mu^ (u = 1/r). (13) 
and collecting terms of equal powers of Muq gives the following two equations: 



#2 



+ (5Mi=3sin^0 (14) 



+ (5m2 = 6(5mi sin (/). (15) 



Solving (14) and (15) for 6ui and 5u2 and substituting them into (12) gives 
the solution 

- = + ^ (l + -Y^) + (iOtt cos - 20</> cos <^ - sm 3<^) . (16) 
Now, one can differentiate (16) and multiply by to obtain 

_ = - _ cos <^ + — sm 20 + cos <^ + - cos 30 + ( - - <^) sm <p) . 

(17) 

After some steps, it follows from (8) and (17) that the bending angle (at = 0) 
to the second-order is given by 

m lf)TT AR^ 

where a = 2'0 of equation (10). The first and second-order mass terms are the 
same as in the expansion in terms of the impact parameter b, see e.g. [13] for 
the Schwarzschild case. 
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2.3 yl-contributions from perturbations of the equations of motions 

Sereno [14] performed calculations based on perturbations of the equations 
of motion in the weak deflection hmit and a term by term integration. He 
considered the light orbital equation from the source at coordinates {rg, ips} to 
the observer at coordinates {ro, 0o = 0}, written in terms of the first integral 
of motion b = ^r^ as 

dr 
r 



1 1 1 2m 

^2 ^2 ^2 y->3 



-1/2 

(19) 



where ta = \/3jA. Next, he expanded the equation for the point of closest 
approach in the weak deflection limit to find 



3rn2 4m3 105m'' 



262 63 2r\ 



(20) 



Then he integrated term- by-term Eq. (19) to find 



4m , / 1 1 \ ISm^TT 128m3 / 1 1 \ SAGBm'^w 
— + 6 - + - ttt; ^ + ^ ^ + ^ 



b \rs ToJ 462 363 5 y^s J 646^ 

3584m-'^ 2m6 m6-'^ / I 1 \ 36'^ / I 1 \ b^ f I 1 



56^ 4 Vr^ ri 40 J Iri 



with the contribution of the cosmological constant appearing in the last term 
and also the term 2bm/r\ that Sereno qualified as local. He then derived for 
small angles the relation between the constant of motion 6 and the angle •& 
between the tangent to the photon trajectory at the observer abd the radial 
direction to the lens as measured in the locally fiat observer's frame 



To 6r3 



^ 6mro 3r, 



62 62r2 



(22) 



He then re- wrote the yl-contribution in a form that, when evaluated at a point 
corresponding to = in Fig. 2 above, reduces to 

yl6=^ 

^^. = -T^> (23) 

in agreement with Rindler and Ishak [5], see Eq. (10) above 

Sereno also proposed in the same paper a derivation based on the lens 
equation and the angular diameter distances where he argues that, to first or- 
der, the lens equation obtained has the same /1-tcrm as one would obtain from 
the integration of the geodesic equations and given by (23). Sereno concludes 
that while the cosmological constant contributes to the bending of light, some 
terms like the term of Ishak and Rindler in Ref. [5] can be incorporated into 
the angular diameter distances. 



Eq. (22) and argues that this term describes the coupling between the lens 



Finally, in another paper [15], Sereno goes back to the term 2bm/r\ of 



8 



and the cosmological constant and since neither the source's or observer's 
position enters into it, this term is clearly local. He argues that contrarily 
to the other yl-term, this local term cannot be integrated into the angular 
diameter distances. 

2.4 The yl-term derived by passing from coordinate angles to physical angles 
and null geodesies integration 

In a first paper [16], Schucker computed the relation between coordinate bend- 
ing angles and physical bending angles and found a result with a A dependence 
in agreement with reference [5] as discussed in section 2.1 above. His calcula- 
tion of the physical angle measured by a static observer on Earth was based on 
using a ratio of proper times of flight of photons from the ratio of coordinate 
times using the Kottler metric. For example, he calculates the physical angle 
Q!T using [16] 



tana7'= — = -^=^= = -^=^= = €t\/ i>(rT) ^ ^Ty ^ — yirx^/S ^ ax 

dr dt^Birr) Ar/-^B{rT) 

(24) 

where er is the coordinate angle. The square root factor term brings in the 

A dependence. 

Schucker did not use any lens equation but used a self-contained method: 
to integrate the null geodesies in Kottler space; to compute area distances to 
the lens and source; and finally to compare his results to lensing cluster SDSS 
J1004+4112, see e.g. [17]. He found, in agreement with [11] that for certain 
clusters of galaxies, the A contribution to bending is not negligible. Schucker 
published a follow up paper on lensing and time delay [18] and another pa- 
per[19] that we discuss in section 3.3 and where he applied his results to an 
Einstein-Strauss (Swiss-Cheese) model. 

3 A light-bending in a Kottler (Schwarzschild-de Sitter) vacuole 

3.1 Kottler Vacuole embedded in an FLRW background 

The cosmological constant has a very negligible effect on small scales but this 
is not the case at large scales of distance between galaxies and clusters of 
galaxies. Therefore, one could measure the term that represent the contribu- 
tion of the cosmological constant to the bending of light using observations 
of large Einstein radii where the lens is a galaxy or a cluster of galaxies, and 
compare them to the first and second order mass terms in the defiection angle. 

The question then is how to model such a system within the cosmologi- 
cal expanding background described by the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) model. As discussed in for example [20], this is still an open 
question since we don't have a stable exact solution to Einstein's equations 
that describes an inhomogeneity (lens) embedded in an FLRW background 
and proposals or alternative approximations have all some sort of caveat. 
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One such proposal is the Einstein-Strauss (Swiss-Cheese) model [23-25] 
with a non-zero A that was used in order to study the A effect on light bending 
in for example [11,19,39] and we use in this section as well. The model is 
known to have an instability to radial perturbations of the boundary of the 
vacuole, see for example [25] and references therein, however, our goal here 
is to find a cut-off location where the Tl-bending of the lens can be regarded 
as accomplished. The model is thus used only to estimate the distance to the 
edge of the Kottler vacuole and use it to give an order-of-magnitude estimate 
of the yl-bending of light term. 

Physically, one is modeling a galaxy or cluster of galaxies and the spacctime 
surrounding it (up to some boundary with an FLRW spacetime) by a Kottler 
vacuole. This Kottler vacuole is in the exact theory embedded into the FLRW 
background using exact matching conditions at the boundary [22,21]. One 
needs to note a point that will have some importance later and that is the size 
of the lens is much smaller than the size of the vacuole which is determined 
from the exact matching between the two spacctimes. In this setting, the rays 
of light enter into the vacuole and propagate inside the Kottler spacetime 
experiencing bending from one boundary of the vacuole to the other and the 
radius at the boundary, r^, enters into the yl-term. 

Mathematically, the Kottler vacuole is exactly embedded into the FLRW 
spacetime using the spacetime matching formalism of [22,21]. The junction 
relations in this case are simple and well-known in the literature [23-25], and 
are given by: 

rt in Kottler = a{t) 1% in FLRW (25) 

and 

Kottler ~ in Kottler ^ Pmatter in FLRW ■ (26) 

Thus, for a given cluster mass, Eq.(26) provides a boundary radius where the 
spacetime transitions from a Kottler spacetime to an FLRW background. One 
shall assume that all the light-bending occurs in the vacuole according to our 
previous formulae, and that once the light transitions out of the vacuole and 
into FLRW spacetime, all bending stops. Noteworthy, the mass-effect falls off 
quickly whereas the yl-effect on the bending of light increases with distance 
from the source; so as we mentioned earlier, the point where to cut off the 
integration, i.e. the boundary of the vacuole, becomes an important factor. 

The results of section 2.1 are applied for the Kottler vacuole here. For the 
small angle (^(, at the boundary, equations (6) and dr/d<p = A give 



1 0b 2m r^/ 20bTO\ 



Inserting (27) into Eq.(8) yields 



277i Ad)hT^ 

ip w tan V' ~ 06 H — 5 ^ + higher-order terms. (28) 

R 6 
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So the bending angle, a, is given, to the smallest order in m/R and A, by 

-«.0-<^,«. — (29) 

Now, to the smallest order, Eq. (27) yields 06 = R/vh, so one can write from 
(29) 

4to ARn 

(30) 

where R is related to the closest approach by Eq. (7) and is the boundary 
radius between vacuole and FLRW, and is determined from Eq. (26). Next, 
this result can be generalized to higher-order mass-terms giving 

m Ihirvn? 305 ARvb 

In [11], the authors used the above setting and some selected systems of distant 
galaxies or clusters of galaxies that are lenses with large Einstein radii and 
found that despite the smallness of A, the Einstein first-order term in the 
bending angle due to these systems is only by some 10^ bigger than the A- 
term and that for the lens systems they considered (see their Table 1), the 
contribution of the cosmological constant term is larger than the second-order 
mass-term, with for example ratios going from 2 to 8 for various cluster and 
galaxy systems. As we discuss in section 3.4 below, Kantowski et al. [39] also 
used the Einstein-Strauss model but found a different yl-term. 



3.2 An upper bound on the value of Cosmological Constant from 
Light-Bending by clusters of galaxies 

While current best constraints on the value of A come from cosmology with 
A ^ 10^^® cm^^, see e.g. [26-36]), it has always been of interest to find upper 
bounds on the value of A from other astrophysical considerations. These in- 
clude upper bounds from planetary precessions of the order of yl < 10~^^ cm~^ 
see for example [37,38] and references therein. The argument used in these 
upper-bound estimates is that the contribution of A cannot exceed the uncer- 
tainty in the targeted measurement. 

As was discussed in [11] and the references therein, the uncertainty in the 
measurements of the bending angle is around Aa ^ 5-10% for several of the 
systems considered in their Table 1. Thus, if the contribution of A cannot 
exceed the uncertainty in the bending angle for these system, then it follows 
that 

'l<|^. (32) 

R n 

For example, with Aa = 10%, one finds from the system Abell 2744 [41,40] 
that 

A < 4.23 lO-^^cm-2. (33) 
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Interestingly, these limits are the best observational upper bound on the value 
of A after cosmological constraints and are only two orders of magnitude away 
from the value determined from cosmological constraints. This result provides 
an improvement of 8 orders of magnitude on previous upper bounds on A that 
were based on planetary or stellar systems, e.g. [37,38]. 

In fact, the limit (33) should be considered as a rough estimate only because 
A also enters into the expression of the angular diameter distances and [11] 
noted that this estimation can be affected by a factor of two or less when 
the uncertainty on A from the angular diameter distances is accounted for. 
Furthermore, there is room for improving this upper bound calculation by 
doing a thorough error analysis where all the other parameters and factors 
that affect lensing along with their uncertainties are measured independently 
from lensing and included in the analysis. 

3.3 Piece-wise integration of null geodesies in an Einstein-Strauss model 

Following on his previous papers on A lensing [16], Schucker performed piece- 
wise integrations of null geodesies [19] in the Einstein-Strauss model with a 
cosmological constant. He calculated deflection angles following his approach 
and Eq. (24) above, taking into account the recession velocity of an observer 
moving with the Hiibble flow of the background spaeetime. As noted in his 
paper, Schucker agreed with the results of [5], that a positive cosmological con- 
stant decreases the bending of light by an isolated spherical mass but found 
that the yl-contribution is attenuated by a homogeneous FLRW mass distri- 
bution around the mass and by the recession of the observer. As he explained 
there, similar screening was also partially obtained in the work of [11] but his 
work using a pieeewise integration of geodesies and various observers allows 
one to quantify this screening. 

Schucker then proceeded with a comparison of the theoretical results and 
lensed quasar system SDSS J1004+4112 [17]. He discussed the dependence of 
the contribution of A to light-bending on estimates of the cluster mass and the 
recession velocities assumed for the observer. Using a Kottler spaeetime, he 
found that a 20% increase of A increases the cluster mass (lens) by 20% for the 
observer at rest and by 10% for the comoving observer (moving with velocity 
V = Hqt). He then redid the calculations in an Einstein-Strauss universe with 
a Kottler vacuole matched inside an FLRW background. He foimd that a 20% 
increase of A decreases the cluster mass by 5% for the comoving observer in 
the FLRW background, noting that the A effect produces more attenuation. 
Schucker concluded that a homogeneous Friedmann background when added 
around the spherical mass and the recession of the observer diminishes the 
effect of A without however canceling it [19]. 

Schucker also points out at the end of his paper [19] that the A bending 
terms can have an interpretation as physical velocity of the observer from the 
lens. He also states there the following: "In their original paper [5] , Rindler and 
Ishak state explicitly that the observer is at rest with respect to the central 
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mass. The same assumption is made in the earhcr hteraturc claiming that the 
bending of hght were independent of A. In my view, this controversy is now 
settled in favor of Rindler and Ishak [5]". 



3.4 Lensing corrections due to A in Flat ylCDM model 

In a recent paper by Kantowski, Chen and Dai [39] the aim is to reexamine 
all higher order corrections to the Einstein bending angle, not just that of A. 
Once again, an exact Swiss-Cheese solution with Kottler in the vacuole and 
flat FLRW around it is used. Angles are measured by comoving Friedmann 
observers. Our interest in this paper is chiefly in what it concludes about the 
yl-correction. Kantowski et al. take careful account of the fact that the radius 
of the Kottler hole in the expanding FLRW universe is larger when the photon 
exits than when it entered. In fact, it is this expansion that in their approach 
gives rise to the most significant part of the correction term. In particular, the 
authors say, without the expansion there is no yl-correction [see after their Eq. 
(23)]. And later, that it arises "since the cosmological constant contributes to 
the extra time ... the Schwarzschild mass has to act on the passing photons." In 
this analysis, the geometric yl-bending inside the vacuole of Rindler and Ishak 
([5,11]) goes undetected. Their lowest order correction in which A appears is 
proportional to 

mVA, (34) 

(their equation (32)) and they find that this could cause as much as a 0.02% 
increase in the deflection angle of the light that passes through a rich cluster. 
Evidently there is need to reconcile these assumptions and findings with those 
of previous authors. But it should also be noted that a term like (??eq:Kterm)) 
can arise simply from aberration in an expanding universe, as we mentioned 
at the end of section 2.1 above. 

Kantowski et al. in their introduction stress the desirability to produce 
formulas which, apart from A, contain only quantities that are measurably 
the same with and without A. "To conclude whether A does or doesn't cause 
bending can easily depend on what is held in common and what property 
is compared in the two experiments" - namely in a universe with and one 
without A. 



4 A lensing and perturbation approaches in an FLRW model 

4.1 Defiection potentials and the cosmological constant 

In Ref. [42], the author shows that the contribution of the cosmological con- 
stant to the light bending angle from [5] can also be derived from gravitational 
potentials, a method that is frequently used in gravitational lensing literature, 
see for example [45,20,44,43]. This method is of course somehow heuristic 
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since it docs not take into consideration some subtleties from the exact the- 
ory and boundary conditions. In this approach, the lens (inhomogeneity) in an 
FLRW background is represented by a Newtonian potential inserted in a post- 
Minkowskian line element or a post-FLRW line element (see for example [45, 
20,44,43]). The metric in such a construction is then given by gab = Vab + hab 
where hab measure the departure from the Minkowskian metric rjab- 

The usual case is that of the linearized Schwarzschild spacetime; hab and 
the associated Newtonian potential <P can be read off the Schwarzschild metric 
in isotropic coordinates as (e.g. [8,12]) 

(1 rn_ \ 2 . 
TTf) + (dr'+r'dn% (35) 

and are given by 

^^_hu^_hu^_rn 

2 2 r ^ ^ 

Next, for a light-ray traveling in the a;-direction, the first-order mass term of 
the Einstein deflection angle is then given by 

Q! = -- / y ±{htt + h^^)dx2 / y±${x,y,z)dx (37) 

J—Xb J —Xb 

where V± = V — V[| is the gradient transverse to the path, see e.g. [45,43]. 
At the vacuole boundary, Xb = y^r"^ — R"^ and the integration gives the usual 
first order result in Schwarzschild 

a..n.,e.„ = ^- (38) 

As we shall see, using this method for the Kottler metric involves some sub- 
tleties. One could argue that in linearized GR and the weak field limit, one 
could add the various hab- It therefore follows that the bending angle can be 
obtained by adding potentials due to different hab- Also, one needs to comply 
with the isotropic FLRW background outside the vacuole so one should also 
use isotropic coordinates for the de Sitter metric, see e.g. [46,47], 

ds^ = ^(L_^] dt^ + L_(^dr^+r^dQ^) (39) 

which is linearized to gab = Vab + hab with, see e.g. [46] , 

htt = -2^=—, hoi=0, hij=-2mij = —dij. (40) 

Applying the same integration of the gradient transverse to the path gives the 
yl-contribution to the bending angle as 



1 /■+^'' f^'"'' ARn 
aA = -^ / ±{htt + h^^)dx = / V^{<^ + ^)dx = — . 



(41) 
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The yl-term (41) found using this method is equal to the yl-term obtained 
from the methods discussed in section 3.1. The method using potentials was 
criticized in [48] which in turn was criticized in [9]. In section 4.2 below it will 
be shown that this term can also be derived from Fermat's Principle. 



4.2 Fermat's Principle and the contribution of A to the lens equation and the 
time delays 

As usual, the deflection angle can also be determined using Fermat's Principle 
and the Euler-Lagrange equations of the variational Principle 

5 J ndl, where n = 1 — 2<?, (42) 

and n is considered as an effective index of refraction of the gravitational field, 
dl is the path of the ray of light. As discussed in for example [45], the null 
curve ds^ = gabdx°'dx'' = for the Schwarzschild potential (36) leads to (42) 
from which the Euler-Lagrange equations give the usual result [45], 

Vi.^dl = (43) 

In [42], Fermat's Principle; was appUcd in order to derive the contribution of 
the cosmological constant to the bending angle and use hab as given by (40). 
Following [42], for a future- directed null curve, 

= = -dt^ (1 - 2^) + (1 - 2^) dl^ = -dt" (1 - + (1 - --) dl^, (44) 



or simply 



dt = {l-{$ + ^))dl = {l + ^)dl. (45) 



It follows from (45) that 



+ + (46) 



Now, following [45], one considers a ray of light traveling along the path dl 
with unit tangent vector e. The deflection angle is given by the change in the 
direction of the null ray. From equations (42) and (46) and the Euler-Lagrange 
equations, it follows that 

_ = - (V(# + ^)-e{e.V{^ + ^)) = -V_l(^ + 0^) (47) 

and 

a A = Bin - Gout = jv±{$ + ^)dl (48) 
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So for a null ray traveling in the x-direction with integration boundary (ri,) as 
given in (41), this approach yields the same yl-term as calculated in sections 
3.1 and 4.1 above, i.e. 

a A = (49) 



4.3 Perturbations and the vacuole model 

Ref. [48] considered the perturbed Robertson- Walker metric where scalar met- 
ric fluctuations are described by scalar potentials, and 

ds'^ = (1 + 2#) dt'^ - a'^it) (1 - 2^) {dx^ + x^dQ'^) 

and where x is the comoving radius, df2 is the line element for the unit sphere, 
and a flat FLRW universe with no anisotropic stresses = <!>) was considered. 

In order to study the deflection inside the vacuole, the authors of [48] aimed 
to find an explicit linear expression for the perturbing potentials responsible 
for the bending of light inside the vacuole so it can be applied in the scheme 
provided in section 4.1. To do that, the authors base their argument on the 
question of how the perturbed FLRW metric compares with the exact Kottler 
metric inside the vacuole. They state that the key to doing this is the transverse 
part of the metric. After a heuristic development, a coordinate transformation 
and some approximations, the authors of [48] arrive at the following expression 
for their function / [in their Eq. (9)]: 



/ = 1 - 2x<P' - 



^ „^ , ainj^l „ 
a in a 



(50) 



The second term inside the square brackets was then dropped out under the 
assumption that the in a perturbed FLRW universe the Kottler vacuole around 
the lens is negligibly small in comparison with the Hubble length and also that, 
in most of the volume of the vacuole the radius r is almost of the same mag- 
nitude as the proper radius of the vacuole. They compare then this expression 
and the Kottler function (2) and find that the potential cf) is given by 

m mr^ 3m 
^=-7-2^+2^ ^''^ 

and has no A term. The authors of [48] conclude that since there is no A terms 
appearing in their potential then there is no yl- terms of the form ylr^/3 in the 
bending angle. They state that while there may be higher order corrections to 
(j) involving A these cannot involve a contribution, and so they claim 

that a contribution to the bending angle from a Ar^ term derived in [5] and 
follow up work is a gauge artifact. The authors state that their solution for 
the potential is a consistant linear solution and the terms neglected were of 
second order. 

On the other hand, Ref. [9] argued that the conclusion reached in [48], 
is the result of too stringent assumptions of smallness that eliminates the A 
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contribution to the bending of light by construction. Ref. [9] argued that it 
was shown in [11,42] that the radius of the vacuole, rf,, is much larger than 
the radius R of the lens itself and when using any expression with the radii 
squared then the corresponding ratio is (^)^ ~ 10'^ or larger, for example, 

for the cluster Abcl2744 [41,40], this ratio is (^)^ w 2500, so they question 
this assumption. They pointed out that this is even more relevant since the 
bending due to the mass is achieved quickly in the vicinity of the lens and 
then changes very slightly as one moves away from the lens to the boundary 
of the vacuole, whereas the effect of A inside the vacuole accumulates up to 
the boundary. This makes the larger size of r& significant. Ref. [9] also argued 
against the assumption of a vacuole size negligible compared to the angular 
diameter distances involved in the lens equations (of Hubble scale) because 
the ratio of the size of the vacuole compared to the Hubble radius is in fact 
about the same order of magnitude (10~^ — 10"'') as the ratio between the A 
contribution term to the first-order Einstein angle. So by neglecting the size of 
vacuole, the effect of A has been neglected with it. Mathematically, [9] showed 
how not discarding the diameter of the vacuole compared to the Hubble radius 
affects the calculation of [48] and restores the contribution of A to the bending 
angle. They did not ignore the terms proportional to <1> in the square brackets 
of Eq. (50) and used the 00-equation of the perturbed Einstein equations (see 
for example Eq. (5.27) in [49]) 

fc2^ + 3-(S- = AnGa'^pmSm, (52) 
a \ a J 

at the same level of approximation where higher order terms on the LHS of 
(52) were kept so the approximation 91n|<?|/91na = —1 used in [48] cannot 
anymore be invoked. [9] modeled a small departure from this relation by —lie 
and followed the same procedure as in [48] finding a potential 

m mr^ Ar^ , . 

plus smaller terms. Their integration gave exactly the 2 mass terms of reference 
[48] plus A terms with a leading term that is precisely — 4f- This A term in 
the potential gives the same A contribution to the light-bending angle (see for 
example [42]) as the one derived from different approaches discussed earlier. 

Other different arguments where made in both [48] and [9] with some 
questions open to debate and the reader is referred to the two papers for more 
discussions. 



4.4 Perturbations of the null geodesies in a first order perturbed McVittie 
metric 

Park used in [50] an approach starting with a perturbed version of an exact 
solution due to McVittie [51] where the solution represents a mass embedded 



17 



in an expanding universe. The metric is given in the notation of [50] by 

fl^Vdf + il + liMtfdX^ (54) 
V 1 "h /X/ 



where X is the comoving coordinate, 



^ = 4a(t)|X-Xo| ' ^^^^ 

with Xq being the loeation of the Schwarzsehild mass and the a{t) the scale 
factor. When the mass vanishes, the metric becomes the FLRW metric and 
when a(t) = 1 the metric reduces to the Schwarzsehild metric in isotropic coor- 
dinates. Park noted that in a real lensing situation in an FLRW background, 
the lens, the source and the observer are all moving according to Hubble's 
law so he proposed to define "physical" spatial coordinates by x = e^*X, 
because the actual distance should be given by (scale factor) x (comoving dis- 
tance) where he used a scale factor a{t) = e^* with H = A/3, in a universe 
driven by cosmological constant (a de Sitter expansion factor). Next, he noted 
that in weak lensing, m in length units is much smaller than any other length 
scale under consideration, and he will work only up to 0(rn), and therefore 
he continued his work in the customized first order perturbed McVittie metric 
given by 

ds^ = - fl V^^-ff^^f)2 ^ 

^ y^{x + e^*q)2 + j/2 _|_ ^2 / \ {x + e^*(j)2 + j/2 ^ ^2 ^ 

(56) 

where he aligned the coordinates to put the lens on the x-axis and notes 
that the cosmic expansion is implemented in the coordinates themselves. The 
Hubble motion of lens with respect to the observer can be seen directly from 
the metric using x = (— e^*g, 0, 0), where g is a constant related to the location 
of the lens at a given time. 

Next, Park wrote the null geodesies equations for this first order metric, 
for the coordinates x and y and integrated them to first order in m and using 
some simplifying assumptions. He used the results to calculate, again to 0{m), 
the angular location of the image of a source S observed by an observer 



, -1 y 

-tan — 



X' 

m 



(57) 

X ys {^^ ^'^sxi. + (l + ( - (1 - Hrs)xi, + ^/(x^^^O^^^"^^^sj^LjM^ ^ 

^Hxj^i^ - Hrs) log 7- T—— r \ + 0{m ) , 

(1 - Hrs)(rs - xsjxj, J 

where /3 = tan~^ || is the undeflected location of the image, and xj, = 

rij/{l — Hrj^). Park used the approximation of a small H and the identification 
of some expressions as being the angular diameter distances in order to re- write 
the lens equation as 
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where 

Park concluded that his result is not in accord with that of [5] where there 
should be a term 0{A) - 0{H'^). 

In [9], the authors questioned the final result of Park since other terms 
including = yl/3 terms were apparently dropped out the calculation at 
some point, leading to the conclusion that A does not contribute to lensing 
except via the angular diameter distances. Ref. [9] claims that even at the 
same level of approximation used in Park's paper, A contribution terms were 
simply omitted from the final lens Eq. (58) including a term with = A/3, 
for reason of smallness. These include terms like 

-/2 



^ I -H'^ + ...]. (60) 



2mdsL a ( Ty2^S<i' 

dsdi 

as communicated by Park [9] , which seems to question the conclusion reached 
by Park. Alternatively, [15] states that the yl-term of [5] is incorporated in the 
angular diameter distances so it can't be found separately in the treatment of 
Park leading to an apparent disagreement. 



5 Other approaches to the A contribution to light bending 

5.1 The Weyl focusing 

Since some early work by R.K. Sachs [53], it is well known that the propa- 
gation of light in spacetime can be rigorously described using the theory of 
geometrical optics in General Relativity [53,54], The distortions of light bun- 
dles are described by the optical scalar equations driven by Ricci and Weyl 
focusing [53,54,45]. The Ricci focusing is given by [53,54,45] 

n = Rabk^k^ (61) 

and the Weyl focusing (up to a phase factor) by 

T = Ca,ibjk''kHH\ (62) 

where Rab is the Ricci tensor, Caibj is the Weyl curvature tensor, and k"' and 
are null vectors and where the notation of [54,45] was used. 
In the Kottler space {Rab = 0), so the contributions of the mass and A 
to the bending of light come from the Weyl focusing as covariantly defined in 
(62). Now, for the Kottler space, the non- vanishing components of the Weyl 
tensor are 

_ 2m 

'^trtr — 

r 

Ctete = -(1 ^) = Ct4,t^/sm^{9) 
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'^rere = 7Z2- = ^r<^r<i>l Sill (t/j 

1 — 2m/r — 

Ce^e^ = 2rm sin^(0). (63) 
Following the work of [54], it is easy to verify for the Kottler metric (1) that 

TZ = and T= (64) 

(equations (27) and (35) in [54]). In Schwarzschild space, h is the impact 
parameter but in Kottler space (not asymptotically flat) h is the constant of 
motion J/E where J and E are respectively the momentum and the energy 
of the photon. Now, we recall that the relation between the point of closest 
approach ro and the constant of the motion h (or b in other notations) is given 
by (e.g [55,38]) 

/i"' ro 3 

For Schwarzschild, the solution ro{h) of this equation is given on page 145 of 
[55]. This solution can be immediately generalized to include A and reads 



ro = — . cos (- arccos ( — 3m\ + A]\ (66) 

.j^:;:^ V3 v v/i^ jj 

This is in agreement with [56] and for Tl = the solution reduces to Eq. 
(6.3.37) in [55] for Schwarzschild. 

Now, if one chooses to express the Weyl focusing fully in terms of the point 
of closest approach ro then (65) and (64) yield 

_ 3m/i^ 3m/ 2m Arn\-^ , , 

So A is present in the final rcsiilt. Similarly, if one chooses to express the 
focusing fully in terms of ft,, then using (66) in order to eliminate ro from (64) 
gives the Weyl focusing, T , as a function of m, /t, and A. In a way or the other 
the final expression for the Weyl focusing has a yl in it. 

This general result is independent of any angle calculations and supports 
independently the conclusion that A does contribute to the bending of light 
in the Kottler space in an invariant and thus coordinate-independent way. 



5.2 Using Friedmann- Robertson- Walker coordinates 

In Ref. [57] , the authors considered an invariant built out of the light frequency 
and the small angle between two incident rays that can both in principle 
be measured at the point of observation. They compared this invariant in 
the de Sitter coordinates and the Friedmann-Robertson- Walker coordinates 
arriving to the conclusion that the A light-bending terms drop out of the 
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lensing equation. However, the authors state in the penuhimate paragraph of 
their paper the following: "Let us note that some corrections on the order of 
~ AvrTo to the lensing effects may exist, as well as other cosmological 
corrections in the general case of the FLRW Universe. However, such "short- 
distance" phenomena are perhaps too small to be of practical interest." 

It was argued in [9] that [57] did an expansion of their Eq. (16) keeping only 
the first term in the invariant considered there and dropping the second term 
of that equation, i.e. wA/O where A = y'yl/S. The neglected term in question is 
ArgTol'^ and represents a contribution of A to the squared invariant that was 
discarded by assumption. 

6 Concluding remarks 

As initiated in the original paper of this series [5] and supported by several 

other papers [14,52,15,16,18,19], although the cosmological constant drops 
out of the coordinate orbital equation of light rays in Kottler spacetime, it does 
nevertheless contribute to the bending of light through the metric itself which 
determines the geometric and measurable properties of the coordinate orbital 
equation. This corrects a long-standing misconception that the cosmological 
constant does not affect the observed light bending angle. Some references 
[14,15,19] suggest that there can be an apparent disagreement between some 
papers due to the interpretion of the A lensing terms. For example [14, 15] 
assert that the A term discussed in the above sections can be viewed either 
as an additional term to the bending angle or can be incorporated into the 
angular diameter distances. [15] also states that this is the reason why the A- 
term was not found separately in [50,57]. Other papers [9] have argued that A- 
terms were missed because of too stringent assumptions of smallness and were 
neglected in the perturbation expansions. In sum, it seems fair to assert that 
the fact that A affects the bending of light has by now been well established, 
while putting the A light-bending terms into a cosmological context is still 
subject to some interpretation and requires further work and clarification. 
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